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Using the effective potential approacfi for composite operators, we fiave formulated a general 
metfiod of calculation of the truly nonperturbative Yang-Mills vacuum energy density (the Bag 
constant apart from the sign, by definition). It is the main dynamical characteristic of the QCD 
ground state. We define it as an integration of the truly nonperturbative effective charge over 
the nonperturbative region (soft momentum region). It is free of all types of the perturbative 
contributions, by construction. For the considered truly nonperturbative effective charge it is finite, 
negative and it has no imaginary part (stable vacuum), as well as it is a manifestly gauge-invariant, 
i.e., not explicitly depending on the gauge-fixing parameter. A nontrivial minimization procedure 
makes it possible to determine the Bag constant as a function either of the mass gap, which is 
^— ■>) ' responsible for the large-scale structure of the true QCD vacuum, or of the effective scale, which 

' separates the nonperturbative region from the perturbative one. We have also argued that the Bag 

, constant is a quant of the energy density which can be released from the QCD vacuum, which, in 

^ ■ its turn, is considering as an infinite and permanent reservoir of energy. 
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I. INTRODUCTION 



(N 

> . 

\l I The Lamb shift and the Casimir effect are probably the two most famous experimental evidences of zero-point 
[ energy fluctuations in the vacuum of Quantum Electrodynamics (QED) [1, 2]. The both effects are rather weak, since 
the QED vacuum is mainly perturbative (PT) by origin, character and magnitude (the corresponding fine structure 
' constant is weak). However, even in this case there have been already made attempts to exploit the Casimir effect in 
, order to release energy from the vacuum (see, for example Refs. [3, |3] and references in the above-mentioned reviews 
' P, Q)- Let us also note that in Ref. [3| by investigating the thermodynamical properties of the quantum vacuum it 
, has been concluded that no energy can be extracted cyclically from the vacuum (see, however Ref. 2] and references 

O • therein). 

I " Contrary to the QED vacuum, the vacuum of Quantum Chromodynamics (QCD) is a very complicated confining 
medium and its dynamical and topological complexity means that its structure can be organized at various levels: 
I I classical [1,0,9 and quantum Q. It is mainly nonperturbative (NP) by origin, character and magnitude, since the 
• • , corresponding fine structure constant is strong. So the idea to exploit the NP QCD vacuum in order to extract energy 
. ^ from it seems to be more attractive. However, before to discuss the ways how to extract, it is necessary to discuss 
K> which minimum/maxmimum amount of energy at all can be released in a single cycle. The vacuum of QCD due 
5^ ■ to its above-mentioned dynamical and topological complexity contains many different components and ingredients. 
[ They contribute to the truly NP vacuum energy density (VED), one of the main characteristics of the QCD ground 
state. It is well known that the VED in general badly diverges in quantum field theory, in particularly QCD. Thus 
the main problem is how to correctly define the truly NP VED which should be finite, negative and it should have 
no imaginary part (stable vacuum). Being one of the main physical characteristics of the QCD vacuum, the truly NP 
VED should not depend explicitly on the gauge-fixing parameter as well. 

Let us also remind that in order to calculate physical observables from first principles in QCD it is not enough to 
know its Lagrangian. It is also necessary and important to know the true structure of its ground state. It is just 
the response of the true QCD vacuum which substantially modifies all the QCD Green's functions from their free 
counterparts. Just these full ("dressed") Green's functions are needed for the above-mentioned calculations. 

The main purpose of this paper is twofold. First of all, to define correctly the truly NP VED, that is free of all 
types of the PT contributions (" PT contaminations" ) , as well as to make it a manifestly gauge-invariant quantity (not 
explicitly depending on the gauge- fixing parameter). Secondly, to discuss in general terms which amount of energy 
can be extracted from the QCD ground state, which in its turn is considering as an infinite and permanent reservoir 
of energy. 
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II. THE VED 



The quantum part of the VED is determined by the effective potential approach for composite operators [lO|, (llj . 
In the absence of external sources the effective potential is nothing but the VED. It is given in the form of the loops 
expansion, where the number of the vacuum loops (consisting in general of the confining quarks and gluons properly 
regularized with the help of ghosts) is equal to the power of the Plank constant, h. 

Here we are going to formulate a general method of numerical calculation of the quantum part of the truly NP 
Yang-Mills (YM) VED in the covariant gauge QCD. The gluon part of the VED which to-leading order (log-loop level 
~ h) is given by the effective potential for composite operators as follows [lOj : 



ViD) ^'-j -^JriHD^^D) {D^'D) + 1}, (2.1) 

where D{q) is the full gluon propagator and Do{q) is its free counterpart (see below). The traces over space-time 
and color group indices are understood. Evidently, the effective potential is normalized to V{Do) = as usual. 
Next-to-leading and higher contributions (two and more vacuum loops) are suppressed by one order of magnitude in 
powers of h at least, and thus are not important from the numerical point of view. 
The two-point Green's function, describing the full gluon propagator, is 



Df^^iq) = - {T^,{q)d(^q\S.) + iL^,,{q)] 1, (2.2) 

where ^ is the gauge-fixing parameter and Tp^i,{q) = g^j/ — q^iqvjq^ = .g/ji> — Lfii,{q). Its free counterpart 
Do = D^uil) is obtained by putting the full gluon form factor d{—q^,£_) in Eq. (2.2) simply to one, i.e., 
= ~ {^i^i^il) + ^L^ti^il)} (l/*?^)- Ill order to evaluate the effective potential (2.1), on account of Eq. (2.2), we 
use the well-known expression, namely 



Tr\n{Df^^D) = 8 x A\ndet{D^^D) = 32 ln[(3/4)d(-q^ + (1/4)]. (2.3) 

It becomes zero (in accordance with the above-mentioned normalization condition) when the full gluon form factor is 
replaced by its free counterpart. 

Going over to four-dimensional Euclidean space in Eq. (2.1), and evaluating some numerical factors, one obtains 
(e, = ViD)) 



dq^ q^ 



ln[l + 3d(g^C)]--d(g^e) + a 



(2.4) 



where constant a = (3/4) — 2 In 2 = —0.6363 and the integration from zero to infinity is assumed. The VED eg derived 
in Eq. (2.4) is already colorless quantity, since it has been already summed over color indices. It does not also depend 
explicitly on the unphysical (longitudinal) part of the full gluon propagator due to the product (D^^D), which in its 
turn comes from the above-mentioned normalization to zero (see above). 

However, it still suffers from the two serious problems. The coefficient of the transversal Lorcntz structure (the 
effective charge d(q'^,£^)) may still depend explicitly on ^. Also, it is badly divergent at least as the fourth power of 
the ultraviolet (UV) cutoff, and therefore suffers from different types of the PT contributions (see below). 



III. THE TRULY NP VED 



In order to define the VED free of all the above-mentioned problems, let us make first the identical transformation 
of the full gluon form factor (which is nothing but the effective charge) in Eq. (2.2) as follows (Euclidean signature 
already) : 



rf(9^0 = diq^O - d^'^iq^O + d'-'iq^O = d^ (q' , + d""' {q\ ^ (3.1) 

where d^'^{q^,£,) correctly describes the PT structure of the effective charge d{q'^,£^), including its behavior in the 
UV limit (asymptotic freedom fl^), otherwise remaining arbitrary. On the other hand, c?^^((7^,^), defined by the 
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above-made subtraction, is assumed to correctly reproduce the NP structure of the effective charge, including its 
asymptotic in the deep infrared (IR) limit, underlying thus the strong intrinsic influence of the IR properties of the 
theory on its NP dynamics. Evidently, both terms are valid in the whole energy /momentum range, i.e, they are not 
asymptotics. Let us also emphasize the principle difference between and d'^ ^ {q^ , £,) . The former is the NP 

quantity "contaminated" by the PT contributions, while the latter one being also NP, nevertheless, is free of them. 
Thus the exact separation between the truly NP effective charge d^^{q'^,£^) and its PT counterpart d^^{q'^,£,) is 
achieved (for the concrete example considered below this separation is unique as well). 

There is also another serious reason for the subtraction (3.1). The problem is that the UV asymptotic of the full 
effective charge due to asymptotic freedom may depend on the gauge-fixing parameter, namely to leading order 



diq',0 -g^^oo inqyAlcDW"^'", (3-2) 

where the exponent Co/&o explicitly depends on the gauge- fixing parameter via the coefficient Cq [l^l, and Aqcd 
the QCD asymptotic scale parameter. Evidently, in the decomposition (3.1) precisely the PT part of the full effective 
charge will be responsible for this explicit dependence on the gauge choice. Subtracting it, we will be guaranteed that 
the remaining part will not explicitly depend on the gauge- fixing parameter (see below). 

Substituting the exact decomposition (3.1) into Eq. (2.4) and doing some trivial rearrangement, one obtains 



dq^ q^ 



ln[l + 3d--(g^e)]-^rf""(9^e) 



(3.3) 



where we introduce the following notation 



epT = J dq^ q^ 

It contains the contribution which is mainly determined by the PT part of the full gluon propagator, d^"^(g^,^). The 
constant a also should be included, since it comes from the normalization of the free PT vacuum to zero. However, this 
is not the whole story yet. The first term in Eq. (3.3), depending only on the truly NP effective charge, nevertheless, 
assumes the integration over the PT region (up to infinity), which should be subtracted as well. The problem is that 
the first term in Eq. (3.3) may still be divergent in the UV limit. A priori, there is only one restriction on its UV 
asymptotic, namely do not contradict to asymptotically free behavior of the full effective charge, to leading order 
shown in Eq. (3.2). For example, if d^^{q^,^) is a linearly falling at infinity function, then the integral in Eq. (3.3) 
still diverges as mentioned above. 

If we separate the NP region from the PT one, by introducing the so-called effective scale q^^^ explicitly, then we 
get 



ln[l 



l + 3d^^(g2,0 



4 



('z^€) 



(3.4) 



dq^ q^ 



Hl + M^^iq\0] 



jNP 



(3.5) 



where evidently 



tprp 



dq' q^ 



ln[l + 3d^^^(g^0] 



iNP 



(3.6) 



This integral represents contribution to the VED which is determined by the truly NP part of the full gluon propagator 
but integrated out over the PT region. Along with epT given in Eq. (3.4) it also represents a type of the PT 
contribution into the gluon part of the VED (3.5). This means that the two remaining terms in Eq. (3.5) should be 
subtracted by introducing the truly NP YM VED evM as follows: 



^YM — £g ^ £PT — e'pT' (3-7) 

where the explicit expression for eyM is given by the integral in Eq. (3.5) (see below as well). Concluding this section, 
let us only note that the above-mentioned necessary subtractions can be done in a more sophisticated way by means 
of ghost degrees of freedom 9||. In other words, the ghost degrees of freedom are to be included into the PT parts of 
the VED. Evidently, the subtracted terms are of no importance for our present consideration. 
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IV. THE BAG CONSTANT 

The Bag constant (the so-called the Bag pressure) is defined as the difference between the PT and the NP VED 
p^ . So in our notations for YM fields, and as it follows from the definition (3.7), it is nothing but the truly NP YM 
VED apart from the sign, i.e.. 



Bym — —iYM — ^PT 



^PT 



1 , 2 2 

/ dq q 







ln[l+3d^^(<z^e)]-|rf^^(<z^c) 



(4.1) 



This is a general expression for any model effective charge in order to calculate the Bag constant (or the truly NP 
YM VED apart from the sign) from first principles. It is our definition of the truly NP YM VED and thus of the Bag 
constant (apart from the sign) as integrated out the truly NP effective charge over the NP region (soft momentum 
region, < (7^ < 9e//)- It is free of all types of the PT contributions, by construction. In this connection, let us 
recall that eg is also NP, but "contaminated" by the PT contributions, which just to be subtracted in order to get 
expression (4.1). 

Comparing expressions (2.4) and (4.1), one comes to the following prescription to derive Eq. (4.1) directly from 
Eq. (2.4). 

(i) . Replacing d{q^) d^P{q^). 

(ii) . Omitting the constant a which normalizes the free PT vacuum to zero. 

(iii) . Introducing the effective scale q^j^ which separates the NP region from the PT one. 

(iv) . Omitting the minus sign for the Bag constant. 

At this stage the Bag constant (4.1) is free of all types of the PT "contaminations" , and it is finite quantity, for sure. 
All other its properties mentioned above (positivity, no imaginary part, etc.) depend on the chosen effective charge. It 
is worth emphasizing once more that the Bag constant and hence the truly NP VED is the main physical characteristic 
of the true QCD ground state. Let us note that in defining correctly the truly NP YM VED (or equivalently the Bag 
constant), the three types of the corresponding subtractions have been introduced. The first one - in Eq. (3.1) at the 
fundamental gluon level and the two others - in Eq. (3.7). For the general discussion of necessity in such types of the 



subtractions see Ref. [14 1 



The quantum part of the total truly NP VED at log-loop level is 

tt=eYM + Nfeq, (4.2) 

where tq is the truly NP quark loop contribution. It is an order of magnitude less than tyM because of much less 
quark degrees of freedom in the vacuum, and it is positive because of overall minus due to quark loop. Evidently, in 
terms of the Bag constant, one obtains 

et = ~BYM[l-Nfvl (4.3) 

where we introduce Eg = vByu and <C 1. 

V. NUMERICAL EVALUATION OF THE BAG CONSTANT 

Eq. (4.1) is the main subject of our consideration in what follows. The only problem remaining to solve is to 
choose such truly NP effective charge d^^{q^,£,) which should not explicitly depend on the gauge-fixing parameter 
^. Here it is worth emphasizing that the implicit gauge dependence is not problem. Such kind of the dependence 
is unavoidable in gauge theories like QCD, since fields themselves are gauge- dependent [T5|. For different truly NP 
effective charges d^^{q'^,£^) we will get different numerical results. That is why the choice for its explicit expression 
should be physically and mathematically well justified. 

For concrete numerical evaluation of the Bag constant (4.1), let us choose the truly NP effective charge as follows: 

d^^iq') ^ ariq') = ^Ip/<i\ (5.1) 

where Kmp is the mass scale parameter (the mass gap) responsible for the large-scale structure of the true QCD 
vacuum. First of all, in continuous QCD it is the NP solution of the Schwinger-Dyson (SD) equation for the full gluon 
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propagator (and references therein). It is well known that in this theory it leads to the linear rising potential 
between heavy quarks, "seen" by lattice QCD 16] as well {{q'^)^^-type behavior for the full gluon propagator). 
Secondly, it does not depend explicitly on the gauge choice. Also, the separation between the truly NP and PT 
effective charge is not only exact but is unique as mentioned above, since the PT effective charge is always regular 
at zero, while the truly NP effective charge (5.1) is singular at origin. Let us note in advance that it leads to many 
other desirable properties for the Bag constant (see below). 

In order to proceed to the explicit numerical calculation of the Bag constant for the truly NP effective charge (5.1), 
let us introduce now dimensionless variable and parameter as follows: 



(z) = 1/z. 



(5.2) 



From the general expression for the Bag constant (4.1) one then gets the dimensionless truly NP YM effective potential 
defined at a fixed effective scale squared g^y^ as follows: 



YM 



>c) 



£YM{qf 



Zc) = -^Z, 



-2 

2 



dz Z 



^d^^(z)-ln[l + 3d''-^(z)] 



NP/ 



(5.3) 



which obviously makes it possible to factorize the scale dependence in the truly NP YM VED (4.1). Performing 
almost trivial integration in this expression, one obtains 



^ym{zc) 



9 In 



(-1) 



z^n 1 



(5.4) 



It is easy to see now that as a function of Zc, the effective potential (5.4) diverges as z^^ at small Zc (the unphysical 
regime, Anp oo because the mass gap is either finite or zero, i.e., it cannot be infinitely large). It converges as 
~ —z~^ at infinity (PT limit, Anp — » 0), i.e. it approaches zero from below. Let us emphasize that in the PT limit 
the truly NP phase vanishes and the PT phase survives only. At a fixed effective scale g^^j it follows that Zc ^ oo, 
indeed. In other words, at a fixed effective scale one recovers correct PT limit, while at a fixed mass gap it evidently 
cannot be recovered. 

From the above one can conclude that the effective potential (5.4) as a function of z^ has a minimum at some 
finite point (local minimum), so its minimization makes sense. In gen eral, by taking first derivatives of the effective 
potential one recovers the corresponding equations of motion [13, UH- Thus the above-mentioned minimization of 
the effective potential (5.4) makes it possible to fix the constant of integration which is Zc in our case. Requiring 
d^YM{zc)/dzc = 0, yields the following "stationary" condition 



z, = 41n[l + (z,/3)]. 



Its numerical solution is 



(5.5) 



mm or) 

Zc ~ ^-^i 

SO at the "stationary" state the effective potential (5.3) can be written down as follows: 



(5.6) 



\Iym{z^^ ) 



1 

2^ 



4^ c / 1 ™„ 



The truly NP YM VED (5.3) and hence the Bag constant thus becomes 



= -0.0263. 



(5.7) 



B 



YM 



^evM = 0.0263qjff = 0.1273 x A 



^NP^ 



where the relation 



(5.8) 



qlff = z—A^p = 2.2A^p (5.9) 

has been already used. 

So, we have explicitly demonstrated that the Bag constant (5.8) at the "stationary" state is finite, positive, and it 
has no imaginary part, indeed. It depends only on the mass gap responsible for the truly NP dynamics in the QCD 
ground state, or equivalently on the effective scale separating the NP region from the PT one. Evidently, only the 
minimization of the effective potential (5.4) makes it possible to fix exactly the numerical relation (5.9) between them. 
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VI. NUMERICAL RESULTS 

In order to complete the numerical calculation of the Bag constant (5.8) all we need now is the concrete value for 
the effective scale Qef / , which separates the NP region from the PT one. Equivalently, the concrete value for a a scale 
at which the NP effects become important, that is the mass gap Ajvp, also allows one to achieve the same goal. 

If the PT region starts conventionally from 1 GeV, then this number is a natural choice for the effective scale. It 
makes it also possible to directly compare our values with the values of many phenomenological parameters calculated 
just at this scale (see Appendix A below). We consider this value as well justified and realistic upper limit for the 
effective scale defined above. Thus, we put 

9e// = ge'// = 1 GeV, (6.1) 
then for the mass gap, on account of the relation (5.8), we get 

Aatp = 0.6756 GeF. (6.2) 

Equivalently, the numerical value of the mass gap Ajvp has been obtained from the experimental value for the pion 
decay constant, F.^ = 93.3 A4eV, by implementing a physically well-motivated scale-setting scheme [l3,[l3]- The pion 
decay constant is good experimental number, since it is directly measured quantity in comparison, for example with 
the quark or the giuon condensates, etc. For the mass gap we have obtained the following numerical result: 

Atvp = 0.5784 GeF, (6.3) 

which for the effective scale yields 

q^ff = 0.857 GeV, qlf f = Q.Tib GeV"^ . (6.4) 

In what follows we will consider this value as a realistic lower limit for the effective scale. One has to conclude that we 
have obtained rather close numerical results for the effective scale and the mass gap, by implementing rather different 
scale-setting schemes. It is worth emphasizing that the effective scale (6.4) quite well covers not only the deep IR 
region but the substantial part of the intermediate one as well. 

For the Bag constant (and hence for the truly NP YM VED) from Eq. (5.8), one obtains 

Bym = -evM = (0.0142 - 0.0263) GeV'^ 

= (1.84 - 3.4) GeV/fm^ = (1.84 - 3.4) x 10^^ GeV/cm^. (6.5) 

This is a huge amount of energy stored in one crn^ of the NP QCD vacuum even in "God-given" units h = c = 1. In 
order to restore the explicit dependence on fi [lO|, the right-hand-side of Eq. (6.5) should be multiplied by it, which 
numerical value in different units is ?i = 1 x 10~^* J s = 6.6 x 10~^^Gey s. Evidently, using the number of different 
conversion factors (see, for example Ref. [lB| or Particle Data Group) the Bag constant can be expressed in different 
systems of units (SI, CGS, etc.). 
Taking into account that 

1 GeV = 0.15 x IQ-^J = 0.48 x 10"" W, (6.6) 
from Eq. (6.5) one finally gets {1 W = IQ-^kW = IQ-^MW = lO^^GT^) 

Bym = (0.88 - 1.63) x 10^^ GW/cm^, (6.7) 

in familiar units of watt (W). This number still indicates a huge amount of energy stored in one cm^ of the NP QCD 
vacuum especially in comparison with the total nuclear energy produced in one year and which is estimated about 
360 GW only. 
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VII. DISCUSSION 



We have already pointed out that the VED is badly divergent as the fourth power of the UV cutoff, so we can write 



e, = -At;y = -X^B, (7.1) 

which can be always done without loosing generality, and where A is the dimensionless UV cutoff. At the same time, 
wc have obtained that the truly NP VED is finite and it is cym = —B. This means that the further analysis can be 
done in terms of the Bag constant. 

Let us imagine now that we can extract the finite portion cym from the vacuum in k different places (different 
"vacuum energy releasing facilities" (VERF)). It can be done by rim times in each place, where m = l,2,3...k. Then 
the remaining (R) VED becomes 



CR = Cg - CE = eg - rimeYM = Cg + B'^ rim, (7.2) 
m=l m=l 

where, evidently, denotes the extracting (E) VED. 

The ideal case (which, however, will never be achieved) is when we could extract the finite portion of the VED an 
infinitely number times and in an infinitely number places, so we should consider the following limits: 



cr = B lim 

(A,fc,n^)— *oo 



(7.3) 



m=l 



Thus, one can conclude that in the ideal case 



ce = B lim 

(A,fc,n„)- 



k 

771—1 



B X A^, 



(7.4) 



diverges quadratically since due to the limits pointed out in Eq. (7.3), we can approximate k ^ X, n„i ~ A, A — > oo. 
At the same time, the initial VED eg diverges quartically, i.e., €g ~ —B x A^, A — »■ oo (see Eq. (7.1)). This means 
that the remaining VED 



en = eg - eE B ]im [-X'^ + X% (7.5) 

A— ^oo 

will be always equal to the initial VED Cg in the limit pointed out above, i.e., in units of B one has 



eR = eg-eE = eg + 0{l/X'^), A ^ oo. (7.6) 

This will be especially so in the real case when the integer numbers k and rim are very big, but finite. The true QCD 

vacuum is an infinite and permanent source of energy. The only problem is how to extract it in the profitable way. 
This is, of course, beyond the scope of this paper, though there are some preliminary ideas on this and related issues. 

It is instructive to show explicitly that the same relation (7.6) takes place between not only the different VED's 
but between different types of energy themselves. Indeed, the different VED's should be multiplied by volume (V) in 
order to get energy. Then from Eqs. (6.5) and (7.4) in units of [(1.84 — 3.4) GeV] it follows that 



Ee = VeE = -(1.84 - 3.4) GeV^^ x lim V n„ ~ -A^ (7.7) 

jm (A,fe,Tlm)— >00 ^^^^ 

since V/ frrt' ~ A'' always when A goes to infinity. 
In the same units, one gets 
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En = Ve„ = -(1.84 - 3.4) GeF-— x A'' - -A^. (7.8) 
So the difference between them which is nothing but the remaining energy Eji becomes 

Er = Eg-EE = Eg + 0(1/A2), a ^ oo, (7.9) 

in complete agreement with the relation (7.6) as it should be. Let us remind that [(1.84 — 3.4) GeV] units are equal 
to [B fm^] units, so in fact this simple analysis is again performed in terms of the Bag constant. 



VIII. CONCLUSIONS 



In summary, we have formulated a general method as to how to numerically calculate the quantum part of the truly 
NP YM VED (the Bag constant, apart from the sign, by definition) in the covariant gauge QCD ground state, using 
the effective potential approach for composite operators. The Bag constant is defined as integrated out the truly NP 
part of the effective charge over the NP region (soft momentum region), Eq. (4.1). At this general stage the Bag 
constant is colorless (color-singlet), finite, and it is free of all types of the PT contributions, by construction. The 
separation of " truly NP versus PT" contributions into the VED is exact and unique because of the subtraction at the 
fundamental level (3.1), as well as due to all other subtractions explicitly shown in Eq. (3.7). 

For the considered truly NP effective charge (5.1) in addition it is negative, and it has no imaginary part (stable 
vacuum). It is also a manifestly gauge-invariant quantity (i.e., does not explicitly depend on the gauge-fixing parameter 
as it is required). The separation of "soft versus hard" gluon momenta is also exact because of the minimization 
procedure. It becomes possible since the effective potential (5.4) as a function of the dimensionless effective scale Zc, 
which separates the NP region form the PT one, has a local minimum. Our numerical value for the Bag constant is 
in a good agreement with other NP quantities such as the gluon condensate. In order to demonstrate this we have 
used the trace anomaly relation, not applying to the weak coupling solution for the P function (see Appendix A). 

Our method can be generalized on the multi-loop contributions in the effective potential, as well as to take into 
account the quark degrees of freedom. These terms, however, will produce numerically the contributions an order of 
magnitude less, at least, in comparison with the leading log-loop level gluon term. What is necessary indeed, is to 
be able to extract the finite part of the truly NP VED in a self-consistent and manifestly gauge- invariant ways. Just 
this is provided by our method which thus can be applied to any QCD vacuum quantum and classical models. It 
may serve as a test of them, providing an exact criterion for the separation "stable versus unstable" vacua. Using our 
method we have already shown that the vacuum of classical dual Abelian Higgs model with string and without string 
contributions is unstable against quantum corrections [H, . 

The Bag constant calculated here is a manifestly gauge- invariant and colorless finite quantity, i.e., it is a physical 
quantity. That is why it makes sense to discuss its releasing from the vacuum. That's the VED in general is badly 
divergent is not a mathematical problem. This refiects an universal reality. Vacuum is everywhere and it always 
exists. Our Universe in general and our real word in particular is only its special type of excitation due to the Big 
Bang. As underlined above, the vacuum is infinite and hence permanent source of energy. The only problem is how 
to release the finite portion (the Bag constant) and whether it will be profitable or not by introducing some type of 
cyclic process. "Perpetuum mobile" does not exist, but "perpetuum source" of energy does exist and it is the QCD 
ground state. Analogously to the Plank constant h which determines a quant of energy emitted by quantum physical 
systems, the Bag constant determines a quant of energy which can be released from the QCD vacuum. 

A financial support from HAS-JINR Scientific Collaboration Fund and Hungarian OTKA-T043455 (P. Levai) is to 
be acknowledged. 



APPENDIX A: THE TRACE ANOMALY RELATION 



The truly NP VED (and hence the Bag constant) is important in its own right as the main characteristic of the 
QCD ground state. Furthermore it assists in calculating such an important phenomenological parameter as the gluon 
condensate, introduced in the QCD sum rules approach to resonance physics [21]. The famous trace anomaly relation 
[2^ 1 in the general case (nonzero current quark masses m^) is 
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where is the trace of the energy-momentum tensor and G"^ being the gluon field strength tensor while as = /Att. 
Sandwiching Eq. (A.l) between vacuum states and on account of the obvious relation (0|O^^|0} = 4ej, one obtains 

46* = (0|^G;:,G^J0) +5]m"^(0|g^g/|0). (A2) 

/ 

Here et is the sum of all possible independent, NP contributions to the VED (the total VED) and (0|gy(7/|0) is the 
chiral quark condensate. From this equation in the case of the pure YM fields (i.e., when the number of quark fields 
is zero Nf = 0), one obtains 



{Q\^G%GlAQ)=AeyM, (A3) 

where, evidently, we saturate the total VED ej by the truly NP YM VED eyM defined in Eq. (4.f ) and calculated in 
Eq. (6.5) at log-loop level, i.e., setting ej = eym + If confinement happens then the [3 function is always in the 
domain of attraction (i.e., always negative) without IR stable fixed point 12]. Therefore it is convenient to introduce 
the general definition for the gluon condensate (i.e., without using the weak coupling limit solution to the (3 function) 
as follows: 



(G^) ^ -{0\^Gl,GlM = -^^YM = ABym. (A4) 



Such defined general gluon condensate will be always positive as it should be. 
The renormalization group equation for the (3 function 



q'^^=l3{a,{q')) (A5) 
after substitution of our solution for the truly NP effective charge (5.1) yields 



/3(a,(g2)) = -a,{q^), (A6) 
so that (3{as)/as = (}(as{q^)) / as{q^) = —1, and Eq. (A. 4) numerically becomes 



(G2) EE _(0|^g;:,G^J0) = {Q\-G%GIM - ^Bym = (0.0568 - 0.1032) GeV\ (A7) 

where we have used the numerical value for the Bag constant (6.5). To the gluon condensate (A. 7) it can be assigned 
a physical meaning indeed as the global (average) vacuum characteristic which measures a density of the NP gluon 
fields configurations in the true QCD vacuum. 

However, it cannot be directly com par ed with the phenomenological values for the standard gluon condensate 
estimated within different approaches [23|. The problem is that it is necessary to remember that any value at the 
scale Kmp (lower bound in the right-hand-side of Eq. (A. 7)) is to be recalculated at the 1 GeV scale. Moreover, both 
values explicitly shown in Eq. (A. 7) should be recalculated at the same parametrization, which is nothing but the 
ratio I3(as)/as- In phenomenology the standard parametrization for the gluon condensate is 



G2 = ( — G^) = (0|^G" G^:,|0) « 0.012 GeV\ (A8) 

TT TT 

which can be changed within a factor of 2 [l9l |. 

Thus in order to achieve the same parametrization the both sides of Eq. (A. 7) should be multiplied by 4(as/7r). For 
the numerical value of the strong fine structure constant we use as — as{mz) — 0.1187 from the Particle Data Group. 
In addition the lower bound should be multiplied by the factor 1 GeV/Ajvp — 1.57, coming form the numerical value 
(6.3). Then the recalculated gluon condensate (A. 7), which is denoted as G2, finally becomes 
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(0.014 - 0.022) GeV^. (A9) 

Just this numerical value for the gluon condensate should be compared with the numerical value coming from the 
phenomenology, see Eq. (A. 8) above. So there is no doubt that all our numerical results are in good agreement with 
other phenomenological estimates [23j . 
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